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Abstract: We describe a general method for calculating the infra-red limit of physi- 
cal quantities in unitary quantum field theories. Using analyticity of Green functions 
in a complex scale parameter, the infra-red limit is expressed as a contour integral 
entirely in the ultra-violet region. The infra-red limit is shown to be the limit of the 
Borel transform of the physical quantity. The method is illustrated by calculating 
the central charge of the perturbed unitary minimal models and the critical expo- 
nents of if 4 theory in three dimensions. We obtain approximate values for the central 
charge which are very close to the exact values using only a one loop perturbative 
calculation. For ip A theory we obtain estimates which are within the errors of other 
more elaborate approaches. 
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Figure 1: Integration contour C in the cut complex plane. 



1. Introduction and summary 

It is important to get a better understanding of the infra-red limit of quantum field 
theory. Conventionally this is studied by extrapolating perturbation theory from 
the ultra-violet using the renormalisation group. The purpose of this paper is to 
show how to augment this approach with information derived from the analyticity 
properties of Green functions. These properties are well-known. For example, the 
Kallen-Lehmann spectral representation of the two-point function of a scalar field 
shows that it has an analytic continuation to the complex momentum plane. We will 
use this to express the infra-red limit as an integral in the ultra-violet region where 
perturbation theory is applicable. For illustration assume that the physical quantity 
F depends on some distance scale s, and that F(s) is analytic in the complex plane 
with the negative axis cut away, then the contour integral 



vanishes, where the contour C is given in figure 1. This means that we can write the 
infra-red limit Fir = hin^^oo F(s) as an integral 



over the infinitesimal circle C and the cut C\. 

We show below that the integral over the cut C\ vanishes for large p, so that the 
infra-red limit Fir then is given as an integral in the UV region over Co where per- 
turbation theory applies. We call this integral I(p) and get that Fir = lim p ^ 00 /(p). 





1 



I(p) is shown to be the Borel transform of F(s) of order k, where k is related to 
the way we analytically continue in the complex scale parameter. The Borel trans- 
form has also been applied to perturbation theory in p6| where the Borel transform 
in the coupling is used. In this case assumptions (justified by details of the dynam- 
ics) have to be made about the poles of the Borel transform. The essential difference 
compared with our approach is that we transform with respect to the scale parame- 
ter. The analytic properties of Green functions, which are needed to define the Borel 
transform, in this case follow directly from general principles of quantum field theory. 

We will use this method on two examples where the infra-red limit of the physical 
quantity F(s) in question is well known, namely the central charge of the unitary 
minimal models perturbed by the relevant field 0(i,3), and the critical exponents of 
<£> 4 theory in three dimensions. 

We also considered the perturbed minimal models in PJ where this method was 
shown to improve upon renormalisation group improved perturbation theory by a 
factor 2. F(s) is given as the two point function of the energy momentum tensor, 
and we will here address some of the details left out of our discussion in [jj. We first 
show the analytic properties of F(s), and we then obtain a bound on the exact non 
perturbative expression for I(p), and this allows us to obtain an approximate value 
for calculating the central charge. 

We employ this approximation method to calculate the central charge for the free 
fermionic and bosonic theories perturbed by a mass term and the unitary minimal 
models perturbed by the relevant operator 0(i,3). The free massive theories can be 
solved exactly and these then demonstrate how I(p) behaves in an exact case. For 
the perturbed minimal models we first calculate F(s) as a renormalisation group 
improvement of a one loop perturbative calculation and show that this expression 
has the correct limit when m — > oo, where m characterises the minimal model. 
We then show how our method improves upon the perturbative renormalisation 
group result. If we use the largest domain of analyticity for F(s) apparent from the 
spectral decomposition there is something of a surprise. The approximation becomes 
dramatically close to the exact value, which is remarkable as it is only based on one 
loop perturbation theory. 

The other example is if 4 theory in three dimensions, where we calculate the 
critical exponents v and rj. ip 4 theory has been shown to be Borel summable of 
order one in the perturbative expansion in the coupling constant. The perturbative 
series will often at most be asymptotic, if for example a theory is ill defined for 
negative couplings the convergence radius will be zero. The question is then if the 
perturbative series captures all of the physics when it is only asymptotic. This 
will depend on the analytic properties of the correlators, and if the theory is Borel 
summable perturbation theory will contain all information about the theory. We use 
the propagator to define functions F u , F v which have the critical exponents v and 77 
as their infra-red limit. All scale dependence is moved into the coupling using the 
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running coupling so that F u (s) and F v (s) are given as asymptotic series in the scale 
parameter s, and we show their analyticity. In this example the actual calculation 
of the physical quantities v and r] is different form the minimal models 
we cannot use the same approximation. The reason is that in this case we use 
higher order perturbative expressions (5,6 and 7 loops) and one cannot find an exact 
expression for the perturbative running coupling (unlike the one loop case). The 
infra-red values for u, i] are therefore obtained from the limiting value of the Borel 
transform, and to calculate this we use a conformal mapping and Pade approximation 
which is done in the appendix. The estimates we obtain for u, 77 are within the errors 
of other more elaborate approximations. 

The main results are that we can use analyticity in the scale parameter to express 
the infra-red physics as an integral in the ultra-violet region, the infra-red limit is 
then obtained taking the limit in the Borel transformed quantity. We show how this 
works in two examples with well known infra-red limits. 

The outline is as follows: In the next section we will discuss analyticity of cor- 
relators, introduce the contour integral giving the infra-red limit and show that the 
contribution from the cut vanishes for large p. We also show that the contour integral 
is the Borel transform of F(s). In section 3 we define F(s) for the central charge, 
show its analyticity and discuss an approximation method for calculating the central 
charge. An exact bound on J(p) is given. In section 4 we use this approximation 
on the free massive bosonic and fermionic theory together with the minimal models 
perturbed by 0(i,3). In section 5 we discuss the y? 4 theory and introduce a function 
giving the critical exponents z/, rj in the infra-red limit, and obtain estimates of v and 
r\ from the limit of the Borel transform. 

2. Defining the contour integral 

The physical quantities of a quantum field theory can be written in terms of Green 
functions which have well known analyticity properties. As an example consider the 
two point function of a scalar field A{p). The Kallen-Lehmann spectral representation 
writes this as a sum over intermediate states 



where dfi 2 c(p 2 ) is the spectral density. Introducing the complex scale factor s G C 



then shows that (A(sp)A(0)) is analytic in the positive half-plane 3?(s) > 0. 

Generally physical quantities will not be analytic in the positive half-plane but 
in some sector S(a) of the complex plane, defined as 




(A(sp)A(0)} = / dp 2 5(/i 2 ) 



1 



s 2 p 2 + fi 2 



S( a ) = {z = re i4) | < r < 00, 



f <0<f}. 



(2.1) 
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Note that ^ 5(a). We want to calculate the infra-red limit Fjr of a physical 
quantity f(x), where Fm = lirn| a .|_ KX j f(x), and we will denote the UV limit Fuv = 
limui^o f( x )- We are considering functions which have well defined limits in 5(a), 
i.e. limisi-j.Qo F(s) = F IR and lim s ^ F(s) = F uv for s G 5(a). 

In the two examples we will consider the physical quantities f(x) are given as 
functionals of a two point correlator resulting in analytic functions F(s). A large 
class of physical quantities will in this way preserve the analytic structure of the 
Green functions in the theory. 

The scale parameter s can of course be introduced into the physical quantity in 
a number of ways. We will define by F(s) the analytical continuation of f(x) where 
s is introduced as F(s) = f(sx)\ x= i so that F(s) has the expansion F(s) = ^F n s n 
around the origin, and F(s) is analytic in 5(a). Choosing another positive power 
f{s 1 x)\ x= i, 7 > 0, gives the same ultra-violet and infra-red limits, but different 
intermediate behaviour. Let us now introduce the scale parameter so that the opening 
of the analytic sector is 2(tt — e') with e' <C 1. F(s) = F(s a ) is analytic in 5 = 
S(2(it — e')) provided that a = 2 ( 7r ° e /) • We will now express the infra-red limit as an 
integral entirely in the ultra-violet region using analyticity. 

Using Cauchy's theorem the analyticity of F(s) implies that the contour integral 

, 1 , I ds —F(s) = -. 1 r (I ds+ [ ds+ [ ds] —F(s) = 0, 
2(vr - e)i Jc s 1 ' 2(vr - e)i \J Co J Cl J C2 J s 

vanishes in the sector 5, where p G M + and the contour C is given in figure 1; e > e' 
so that C C 5. In the limit where — > oo, the contribution from the contour C*2 
(see figure 1) becomes 

lim 1 / H d9e pe -' B I^F{r IR e w ) = - lim F{s) = -F IR . (2.2) 

rm^oo 2[TT - e) J 7T _ e |s|— oo 

The angular integral and the limit rj R — > oo can be interchanged as the integrand 
in (|2.2|) is bounded by a constant in the limit ri R — > oo. We then get the integral 
representation of Fir 

1 ft e p/s t e p/s \ 

Analogously by considering 

, 1 : / ds —F(s) = , 1 : ( / ds + / ds + / ds) —F(s) = 0, 
2(vr - e)i Jc s 1 } 2(vr - e)i \J Co J Cl J C2 ) s 

we get that 

F ^ - why, H " s t f ^ + 1 ds T F w) • < 2 ' 4 > 
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We denote the integral along the contour C\ (close to the cut) by 



1 f e p/s 
cut(p) = — / ds F(s). (2.5) 



2(n-e)i 



The integrand in this integral is damped by the factor e~ p ^ s \ and we will show that 
linip^oo cut(p) = by showing that it is bounded by a finite integral for all p allowing 
us to take the limit p — > oo in the integrand. We substitute s = re ±l ^~ e ' for points 
on C\ in the upper and lower half-plane, the integral becomes 

— 1 / r uv pPe'^-^/r rr IR pe*(*- e )/r 

cut(p) = — / dr F(re l{n - t] ) + / dr F(re-^- e) ) 

2(7r-e> \ J rtR r J ruv r 

I f riR , e 



7r — e 



/ dr Im[e ipsin( "" e)/r F(re~ i( "- £) )] 

•Jr uv r 
r r iR e pcos(7r— e)/r 

/ dr (cos(p sin(7r - e) /r)Im[F(re" l(7r_e ^ 

J mv r 



+ Re[F(re- iin - e) )} sin(psin(7r - e))) . 

We divide the r interval into (r uv , 1) and (1, r IR ) and write cut(p) = cut uv + cutjR. 
The function F(s) is finite in S hence there exists a constant q > so that in the 
limit ruv ~^ 

|cu%y| < q dr > for p — > oo. (2.6) 



We know that F(s) — > Fir G M for \s\ — > oo from S, hence Im[i ?1 (re - '*^ 7r_e ')] — > 
for r — > oo. If Im[F] falls off like r~ s for some 5 > then in the limit rj£ — ► oo, 
3 fci, k 2l ks > 0: 

/•°° (sm(k 2 /r) + k 3 r~ s ) 
\cut IR \ < h / dr± — \^LJ. i L < oo f or ail p e r 2 .7 

For a general F(s) where the fall off might be slower we keep a finite rjjj then \cutm\ 
is again finite for all p and a finite rm introduces a 0(— ) term in ( |2.3| ) which is 
negligible for rm large. As e < 1 we can replace 2 ( 7T 1 -e)i w ^ n 2^7 §i ven a term 0(e) 
on the right hand side in (|2.3|) which is again negligible for small e. We can then 
define 



1 f e p/s 

I m (p) = lim — / da F(s) (2. 

- 'Co * 



so that Fir — h m p^oo Iir(p) an d we have succeeded in writing the infra-red limit as 
the limiting value of a contour integral in the ultra-violet region. 

We will now use Cauchy's theorem again to rewrite Im{p) and this will show 
that Iir(p) is the Borel transform of F(s). Instead of integrating over Co we will 
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integrate over C given as the path from the origin along the ray arg(s) = —it + e and 
then anti-clockwise along \s\ = x > until arg(s) = 71 — e and then back to the origin. 
It follows that this contour integral is independent of the choice of x and e as long 
as e' < e < ~. The upper limit ensures that the rays stay in the negative half-plane 
(where e P//s is a damping factor), and the lower limit that F(s) is analytic on the 
contour. Using this contour the ultra-violet limit ruy — > can be taken explicitly 
by extending the rays to the origin. 

Using that F(s) = F(s a ) for a = 2 ( 7T a L e ') then amounts to 

where k = -, p = p l ' k and C' is the contour where the rays satisfy |arg(s)| = 
|(tt - e")/k\ for any e" G (e', f ). We will write e" = (tt - e)/2 for any e G (0, vr - 2e'), 
then I arg((^) fc )| > n/2 and the integrand in (|2.9f) is again damped on the rays. Also, 
on the rays is | arg(s)| = (it + e)/2k = f < Oi/2 so that F(s) is analytic on the 
contour, C C S'(a), and the contour integral is therefore well defined. 

Equation (|2.9| ) then shows that Iir(p) is the Borel transform of F(s) of order k 
[ [16|j . The Borel transform of order k of a formal power series h(z) = Y2 n h n z n is given 
by B k (h)(p) = J2 n r(i+n/fc) ' using the inte g ral representation = ^- f & ^ e l/s s n 

it then follows that Q2.9Q is the Borel transform of order of i^(s) by inserting a series 
expansion for -F(s). 

The contour integral in ( |2.9| ) is independent of the contour C', i.e. in the choice 
of e for e G (0, n — 2e'). This shows that equation ( |2.9| ) holds for all k with F(s) = 
F(s 1,/A: ) and a>it/k 1 because we can always find an e > so that a > (it + e)/k, 
which again implies that C C S(a) and | arg(p/s) fc | > tt/2 which again makes the 
contour integral well defined. 

We have then shown that: the infra-red limit of a physical quantity, Fir, is the 
limiting value of the Borel transform of F(s), and changing the way in which the 
scale parameter is introduced amounts to changing the order of the Borel transform. 
The order k has to satisfy the bound a > it/k where a determines the analytic sector 
of the physical quantity F(s). 

One way of calculating the contour integral in ( |2.9| ) is to insert a series expansion 
of F(s) around the origin, but because F(s) is only analytic in a sector this series 
can only be an asymptotic expansion of F(s) 1 . An analytic function has a unique 
asymptotic expansion, but an asymptotic series might be the asymptotic expansion 
of several analytic functions. However in the case where a function is analytic in a 

x If f(z) is an analytic function in a neighbourhood of the origin it follows that its asymptotic 
series around the origin is convergent and equal to the power series of f(z); the lack of analyticity 
at the origin forces the asymptotic series to diverge. A formal power series f(z) — fnZ n is 
asymptotic of order k if 3 C, K > : |/ n | < CK n T(l + n/k), this shows that the Borel transform 
of order k, of such a series, has a non zero convergence radius. 
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sector of opening a > |, which is the situation we have above, the map between the 
analytic function and the asymptotic series of order k is injective (but not surjective) 
[16|j . This means that all information about the exact function is contained in the 



asymptotic series. 

An asymptotic series of F(s) in the scale parameter can be obtained doing the 
perturbative expansion in the coupling and then inserting the running coupling con- 
stant. Assume that F is given as a formal perturbative series, which we write as 
F(s,g) = J2n-Fn(s)g n . The Callan-Symanzik equation states that a theory is in- 
variant, i.e. the correlators are invariant, under a scale transformation if the cou- 
plings change according to the renormalisation group. A scaled quantum field the- 
ory can therefore equivalently be described by a theory on the same scale, but with 
couplings changed according to the renormalisation group. In this way scale de- 
pendence of a theory can be moved into the running coupling g(s). Moving all 
scale dependence into the running coupling we get an asymptotic series in the scale 
F(s = l,g(s)) = ZF n s n 2 . 

The examples we consider below have a = tt — e', where e' < 1, hence from the 
constraint a > k = 1 + 5 for some 5 > 0. We will generally try to minimise the 
order k (k — > 1), thus maximising the analytic sector and the convergence of the 
Borel transform. 

3. Defining F(s) for the central charge 

As the first example we want to calculate the infra-red central charge of a 2 dimen- 
sional unitary quantum field theory. The renormalisation group flow of the central 
charge is governed by Zamolodchikov's c-theorem M. This states that for a uni- 
tary and renormalisable quantum field theory in 2 dimensions there exists a function 
which is monotonically decreasing along the renormalisation group flow, and which is 
stationary only for conformally invariant theories where it takes the value of the Vi- 
rasoro central charge. The c-theorem implies that the infra-red limit, where the scale 
goes to infinity, and the ultra-violet limit, where the scale vanishes, are fixed points 
of the renormalisation group. In 2 dimensions scale invariance implies conformal 
invariance so in these scaling limits we have a conformal field theory characterised 
by the central charge. We define the energy-momentum tensor as in |2| 

where V = Vo^S 1 ™ -1 ) = 2tt in 2D and W[g] is the effective action. 

2 If F(s) is asymptotic of order k! (F(s) <E C[[s]]i/k') we will choose k = k', and then the Borel 
transform of F has a non-zero convergence radius and the analytically continued value at infinity 
uniquely determines Fir. The analytical continuation can be done by a conformal mapping, we 
will discuss this in the appendix. 
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The correlator (T zz (z)T zz (0)) gives the ultra-violet and infra-red central charges 
in the limits z —> and z — > oo respectively. In (T Atiy (x)T p(T (0)) was written using 
the Kallen-Lehmann spectral representation: 

K fj 2~/ 2\ [ d2 P ipx(9fiuP 2 -PvPv){g P oP 2 -PpPa) 



(T^(x)T p(J (0)) = — — / dp c(p 



-e 



3 • 16 J ^ ^ ' J (2vr) 2 p 2 + p 2 

and it follows that 

(T zz (z,z)T zz (0,0)) = — \ dp 2 c(p 2 ) / -JL——p\ (3.2) 



3 -16 J ^ ^ J J (2tt) 2 pp + p 2 

where we use the usual complex variables 3 z, z, and c(p 2 )dp 2 is the spectral density 
which represents the density in degrees of freedom of the quantum field theory at the 
mass p. If we scale z, z by a positive real dimensionless parameter s (s > 0) we get 

7r f°° r d 2 a (4(^+9*0 a 4 

(TUs.MOM - 5^ jf / p^-^fl- (3-3) 

In the ultra-violet limit where s — > then _ , g 2 2 ~~ > ^ an d (13.31) becomes 

qq+s^p^ qq v l k 

(T zz (sz,sz)T zz (0,0))^^J o dp 2 d(p 2 ) = ^ fors^O. (3.4) 

The ultra-violet central charge is therefore C[/y = J °° dp, c(p) . To calculate the infra- 
red limit we first note that 



d 2 q e^ qS+ ^ _ 4 o4 / d 



1 



,0 v> - > 2$ = 2 br g (^> s a0> ( 3 - 5 ) 

(271) qq + s 2 p 2 \oz J 

where G(z,z,p) is equal to the free Bose propagator at mass p. G(z,z,p) can be 
written in terms of a modified Bessel function || 

f d 2 p e^> K (\x\p) 

G{x,p) = = — . 3.6 

J {2tt) p + p 2tt 



Performing the differentiation we can then write (|3.3| ) as 

1 r°° 

(T zz {sz,sz)T zz {0,0)) = 2 48g4z4 J dp 2 c(p 2 )ps\z\((p 3 s 3 \z\ 3 

+ 2Aps\z\)K (ps\z\) + (8Ai 2 s 2 |^| 2 + A8)K l (ps\z\)). (3.7) 

In the infra-red limit where s — > 00 Ko(p\z\s) and K\(p\z\s) have the asymptotic 
behaviour |J e~ MS ' z and the only contribution to ( |3.7| ) comes from the massless limit 
where p — > 0, hence 

(T zz (sz, sz)T zz (0, 0)) -> lim — ^-j jf d/i 2 g(/i 2 ) = for s -> 00, (3.8) 



3 z = + is, z = tE — ix and <i 2 z = e? 2 a; = A di_E = — 5^2 A g?z. 
4 The Fourier transform of (ir/24)p 4 /pp is 1/z 4 . 
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so that cm = lim e ^ J Q dp 2 c(p 2 ). This shows that cuv > cjr as the spectral density 
is positive for a unitary theory, and this is another way of showing the c-theorem 
0. This representation of the central charge using the spectral representation also 
shows that the central charge measures the number of (massless) degrees of freedom 
of the CFT. We define the function 5 

F{s) = 2z\s 4 (T zz (sz, sz)T zz (0, 0)> (3.9) 

z=z=l 

here s and F(s) are dimensionless and s G M+. This function then satisfies 

Hs) -> ( CUV !° r ° + ' (3.10) 
[ cm tor s —>■ oo. 

3.1 Analyticity of F(s) 

We will show that F(s) is an analytic continuation of F(x), x G M+, for s G S = 
S(n — e') (with e' <C 1); to show this write F(s) = Jj^duf^Sju). For F to be 
holomorphic in 5* then f(s, v) must be holomorphic in S for all v G M, and both / 
and ^ must be integrable over the set M. Using ( ^.7| ) above we can write F(s) as 

i r°° 

F(s) = — dfi 2 c{fi 2 )fis ( {fi 3 s 3 + 24 fis)K {^s) + {8 fi 2 s 2 + 48) K^s)) . (3.11) 

K v (z) is holomorphic in S so f{s,v) is clearly holomorphic in S. The modified 
Bessel functions also satisfy that l-fT^z)! is bounded for \z\ > e for any e G M + 
and | arg2;| < |, which is the case for z = [is when s 6 S. For large values of \s\ 
integrability is ensured by J °° dfi 2 c(/x 2 ) = cuv < oo and the asymptotic behaviour 
K v {z) = \p^e~ z (X + 0{z~ x y). Around the origin K (z)z — > for z — > and 



K\{z)z — > 1 for z — * 0, hence -F(s) is integrable. dF } s ' is shown to be integrable in a 
similar way and the analyticity is shown. 

From the form of F(s) in (|3.11|) it follows that F(s) has a limit value F(s) — > <i 



for s — >• from S 1 . From ( |3. 1Q|) it follows that d = Cuv, an d we define F(0) = cuv- It 
also follows from ( |3.11|) that F(s) has a limit value for \s\ — > oo from S 1 : -F(s) — > tf 
and ( P-10| ) again sets d' = cir. 



3.2 The approximation for cjr 

We have the following representation of cir and Cuv 

op/s r p p/s 



\fr ppl s r e p/ s \ 

ciR= — [ / ds F{s) + / ds F(s) = I IR (p) + cut(p) (3.12) 

Zm \JCo s JCi s / 

-I / r e P s r e P s \ 

c uv = 7T-\ / ^— F(s) + / ds—F(s) = I uv (p) + cut(p), (3.13) 



2th V./c 2 s J Cl 



5 When we set \z\ — 1 then fi becomes dimensionless in ([T?]), we also denote this dimensionless 
quantity by fi. 



9 



where in both cases the contribution from the cut is rapidly decreasing in p. Note 
that in these relations the integral is performed in the opposite scaling limit of the 
quantity we calculate. We will here concentrate on ( |3.12|) as perturbation theory can 
be applied to F(s) when s G C , we will call Iir(p) = I{p)- 

We can choose coordinates in the coupling constant space so that the ultra- 
violet fixed point corresponds to g(s) = 0. In the ultra-violet limit where s->0we 
may describe F(s) by perturbation theory as g(s) — > 0. The nth order perturbative 
approximation of F(s) is denoted by F n (s) and the corresponding integral by I n {p)- 
In the limit of large p the contribution from the cut vanishes and we get 

c m « lim J n (p) = lim F n {s) (3.14) 

p^oo s— >oo 



where the last equality follows setting s' = s/p in (|3.12| ) and then taking the limit 
p — > oo in the integrand valid for all Tuy > 0. Moving all scale dependence into 
the running coupling g(s) we can write F n (s) = § n (g(s)). Let g* IR denote the first 
non trivial zero of the perturbative j3— function, i.e. g(s) — > g* IR for s — > oo. Then 
equation Q3.14|) becomes 



lim F n {s) = lim $ n (g(s)) = $ n (g* IR ) = c* IR (3.15) 



-OO 



which is the perturbative estimate of cir we want to improve. 

F n (s) is the ultra-violet perturbative approximation to F(s) and the integra- 
tion range in J(p) is compact so lim p ^ I n (p) = \im p ^ I(p) = cuv- I n (p) therefore 
provides a good approximation to I(p), for small enough p, since the power series ex- 
pansion of I(p) is controlled by the small s expansion of F(s) for which perturbation 
theory applies. This is illustrated in ( |4.30| ) below for the minimal models. For larger 



values of p higher order terms in the expansion of I n {p) become important and the 
coefficients of the expansion of I n {p) and F n (p) part company. If cir < c* IR and if the 
region where I n (p) is a good approximation to I(p) is large enough, then I n {p) will 
have a minimum before approaching its limiting value of c* IR . Since this minimum 
occurs at the largest value of p for which I n (p) is a reasonable approximation to I(p) 
and the true value of cir is given by /(oo), it is this minimum of I n (p) that we will 
use to provide a better estimate of cir. The approximation then becomes 

cir = I n (Pm) (3.16) 

where p m is the value where I n (p) attains its minimum. Below we consider the quan- 
tity Ac exact = cuv — cir, and denote the approximation to it Ac approx = cuv — In{Pm), 
we call the perturbative value Ac pert = cuv — lim p ^oo In{p) = cuv — c* IR . This ap- 
proximation rests upon the assumption that the exact function is monotonically 
decreasing from c^y to cir, or at least that its minimum value is close to cir. Be- 
low we will show that this is indeed the case for k — 4, at least to a very good 
approximation. 
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3.3 Exact bound on Ac 



As discussed above can the exact value I(p) possibly be smaller than the asymptotic 
value Cir. We will denote by p' m the value where I(p) attains its minimum and 
Ac est = cuv — I (p'm) j then Ac exact — Ac est measures the undershoot of the exact 
function I(p) compared with its asymptotic value cir (see figure 3). From the spectral 



representation of F(s) in ( 3.11 ) we can obtain a rigorous lower bound on this exact 



undershoot. In [J2J it was shown that the spectral density can be written as 

c(p 2 ) = c IR S(p 2 ) + c(/i 2 ), 
we showed in (|3.4j ) that f dp 2 c(p 2 ) = cuv hence J dp 2 c(p 2 ) = Ac exact . Using this in 



( |3.11|) we get 

Ac est = c uv - I(p' m ) = c uv - ^— \ 

2m Jc 

= Ac exact - / dp 2 c(p 2 )T(p 2 ,p' m ) 



1 / dS e P'm/s F ^ 



where 



np 2 ,p'J = ^~ I -e p '- /s ^(((ps^f 



2ni Jc s 48 

+2Aps 1/4 )K (ps 1/4 ) + (8(ps 1/4 ) 2 + A8)K 1 (ps 1/4 )) 

and F(s) = ^(s 1 ^) with k = 4. Rescaling s allows us to move all p dependence into 
P'm(p) = P'mP 4 i so that we can write T(p 2 ,p' m ) = T(p' m (p)). Unitarity ensures that 
c(p 2 ) > hence 

POO 

Ac exact - Ac est > mmT(p m (p)) / dp 2 c(p 2 ) } (3.17) 
^ Jo 

Ac est > so that 

^ A ° ex f ~ A ° est > minT(p). (3.18) 

Ac exact p>o 

The lower bound in the relative undershoot of I(p' m ) therefore equals min p T(p). In 
figure |2| we plot T(p) for k = 4. With this value we get that min p T(p) = —0.0232 
so the relative overshoot in Ac est compared with Ac exact is maximally 2.3%. It 
follows from (|3.17|) that the bound in (|3.18|) is only saturated for free theories (where 
c(p 2 ) oc S(p 2 — m 2 )) for general interacting theories will the relative overshoot be 
smaller than | min p T(p)|. We show in figure |^ the type of behaviour that we expect 
and which is confirmed for the minimal models below. 

The actual choice of k = 4 then is a compromise between maximising the range 
of analyticity (small k) and minimising the undershoot in I(p' m ) (large k). 
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Figure 2: Numerical integration of T(p) given in (|3.3|). 



Ac est 




Figure 3: The expected behaviour of cjjv — In(p) & n d cuv ~ I(p)- 



4. Application to different 2D models 



In this section we calculate the central charge for the free bosonic and fermionic the- 
ory perturbed by a mass term. These theories are gaussian and we can calculate F(s) 
exactly, the central charge can then be written as lim p _ +00 /(p) as the contribution 
from the cut then vanishes. The infra-red limit of these theories is trivial (c/# = 0) 
as all the degrees of freedom are massive, and therefore decouple when approaching 
the infra-red fixed point where the scale goes to infinity. The free theories are none 
the less important to consider as we can here obtain the exact function I(p). 

In 4.3 we consider the unitary minimal models perturbed by the relevant oper- 
ator 0(i,3), this theory has a non-trivial infra-red fixed point. We obtain the renor- 
malisation group improved perturbative calculation of cm and compare with the 
approximation ( |3.16| ). 
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4.1 The free boson 

We take the action for the free bosonic theory in 2 dimensions with a mass m to be 

S = J d 2 x (±d lt <p(x)& t <p(x) + ^mV 2 wj • (4.1) 

The perturbation away from the conformal field theory is thus given by the purely 
massive term ^m 2 (p 2 . The theory is still a free theory off criticality and the correlator 
(TT) can be calculated exactly in the whole scaling region from the ultra-violet to 
the infra-red. With this normalisation the energy-momentum tensor becomes 



TO, z) = T zz {z, z) = -2tt : d<p(z, z)d^{z, z) : (4.2) 

with the correlator 

(T(z, z)T(w,w)) = (27r) 2 (: d<p(z, z)d<p(z, z) :: d(p(w,w)d(p(w,w) :) 

= 2(2ir) 2 (d(p(z,z)d(p(w,w)} 2 , (4.3) 



as only the double contractions survive. Using the form ( |3.6| ) of the free propagator 
then (|4.3|) is 2(d z d z K (m\z\)) 2 where we have set w = using translation invariance. 
We now use the identities K' n = — |(iT n _i + K n+1 ), K 2 (x) = K (x) + ~Ki(x), 
K-i(x) = KAx) and ([4.31) becomes 



m 2 \z\ 2 



(T{z,z)T '(0,0)) = ' 4 ' (AK^mlzD+m^z^K^mlzD+AmlzlK^mlzDKoimlzD) , 
hence F(s) with k = 4 becomes 

i^( m / /4 )+»^oK 1/4 ) + ^oK 1/4 )^iK 1/4 )) • 



l/2\ 2 

ms 1 



F(s) 



Knowing F(s) exactly the central charges follow from ( ^.1U| ) directly, but let us use 
the relations ( |3 . 1 3| ) , ( |3.12 ) in the limit of large p where the cut vanishes. For cjjv we 



substitute s' = ms 1//4 , p' = in the contour From ( p.!3| ) we then get that (for 
rm = ^t) 

7T 

I uv (p) = i/' d9e m e"' e216 (Ae-™Kl{e w ) + K 2 (e w ) + Ae~ w K^K^)) . 

47T J_2L 

(4.4) 

This is not expressible in terms of elementary functions, but it can be calculated 
using the analytical properties of the Bessel functions. The integration contour 
can be collapsed into a contour running along the imaginary axis together with 



13 



0.9975 - 



Figure 4: a) Numerical integration of Ijjv(p) given by (4.4) compared with the exact 
value cjjv = 1- b) Numerical integration of Ijr(p). 



an infinitesimal semi-circle around the origin. Setting ix 
contributions 



e t0 we then get the 3 



1 



dx+ / dx+ / dx 

1 Je JC e 



e~ px2 {-AxK^ixf + x 3 K {tx) 2 - Aix 1 ] K {ix)K 1 {ix)) . 



For p — > oo only the contribution from the infinitesimal semi-circle C e survive and 
here we can insert the asymptotic form for Ki(z) ~ z" 1 and K (z) ~ — log 2:. In 
this limit only the first term with K\ will contribute as is seen from the asymptotic 
form. Taking into account that we only integrate over half a circle we then get the 
well known result cjjy = 1, which is an exact result as the contribution from the cut 
vanishes in this limit 6 . To calculate the infra-red central charge we use (|3.14|) hence 
Cir = linip^oo Iir(p) = lim^oo From the asymptotic form of the modified 

Bessel functions K n (x) ~ e~ x , together with ( |4.1| ), it follows that cir = 0. Figure £|a 
shows the exact function Ijjy(p) computed by a numerical integration of ( [4.4] ) using 
a NAG Fortran Library integration routine. The figure shows how the contribution 
from the cut vanishes in the limit of large p. The exact infra-red function Im{p) is 
plotted in figure |4]b, it is important to note that the minimum value of Im{p) is very 
small, namely —0.019. 



°We could also directly have used the analogous of ( p.141) namely linip^oo Iuv(p) = lini s ^ F( s ) 
which is again seen changing variable in (3.13) and taking the limit p — > 00 in the integrand valid 
for all rm < 00. 
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4.2 The free fermion 



The free massive fermion has the action 

S = J d 2 x {^d^ + i)di) + im^) . (4.5) 

A calculation analogous to the Bose case gives, for k = 4, that 

F(s) = ^ (Kl{msV±) (l + -i^) - Kl{m S y^ . (4.6) 

The same arguments as above yields cuv — \, an d the infra-red contribution Cm = 
lim s ^ 00 F(s) = 0, again by using the asymptotic form of the Bessel functions. Again 
we can obtain the exact functions I(p) by a numerical computation. We get the same 
behaviour of Iir(p) as in the bosonic case, but with the minimum value —0.0048, so 
this exact function does again not differ much from the infra-red central charge at 
its minimum value. 

4.3 The unitary minimal models 

The off critical quantum field theory picks out a specific renormalisation group flow 
from the ultra-violet to the infra-red conformal field theory. If the quantum field 
theory is in the neighbourhood of one of the renormalisation group fixed points A* 
in coupling constant space (we may be choose the coordinates so A* = corresponds 
to the ultra-violet CFT) the action can be written as 

N 

S = S CFT + J2 Xl d 2 x$i(x). (4.7) 

i=l J 

Here the $j(x)'s are scaling fields with scaling dimension Aj, the coupling constants 
X 1 then have mass dimension [A 1 ] = 2 — Aj = y^. From (|4.7|) it follows that %ji is 
the renormalisation group eigenvalue of the scaling fields $j. For renormalisable 
quantum field theories we need yi > and the scaling operators therefore have to be 
relevant (y > 0) or marginal (y = 0). For relevant operators we will move away from 
the fixed point when the scale increases and Scft thus corresponds to an ultra-violet 
critical point. 

We consider quantum field theories which have the unitary minimal models 
Ai(m) as their scaling limits. The central charge for Ai(m) is given by 

c(m) = l- , 6 , - , m = 3,4,... (4.8) 
m[m + 1) 

There are m(jn — l)/2 primary fields 

4>(p,q) = 4>(m-p,m-q+l), 1 < P < 171 - 1, 1 < q < TU, 
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with the conformal dimensions 

((m + l)p — mq) 2 — 1 



h M ~ h M 4m(m + 1} , (4-9) 



so the primary fields are scalars (spin zero). It follows from (|4.9|) that 2{m — 2) of 
the primary fields satisfy A = 2h < 2 and are thus relevant operators 7 . We will 
consider the minimal models perturbed by the relevant operator 0(i,3) given by the 
action 



S = S M (m) - X j d X (p(i, 3 )(x). (4.10) 
The reasons for choosing this operator are that: 8 

i) 0(1,3) is a relevant field, ^(1,3) = 1 — < 1, which exists in all Ai(m). 

ii) <P(rn,n) form an algebra under the operator product expansion, from the fusion 
rules of minimal models it follows that the set 4>h, n ) constitutes a sub-algebra in 
which only 0(i,i), <f>n,2) an d 0(i,3) are relevant as seen from ( (191) . 0(i, 3 ) is normalised 
so that the structure constant C(i j 3)(i )3 )(i ) i) = 1 where I = 0(i,i) = 4>( m -i,m) is the 
identity, and it has the self coupling Cn^n^n^ = b{m). It does not couple to 0(i,2) 
so it has no coupling to other relevant operators in the sub-algebra <fin n ) B ||) 0]- 
This means that there is a renormalisation group flow connecting the ultra-violet and 
infra-red fixed points along the direction of 0(i,3) so it is consistent to include only 
the one relevant field 0(i,3) 9 , i.e. it is a geodesic renormalisation group trajectory [^]. 
Hi) 0(i,3) is the least relevant field, and the perturbation in ( |4.10| ) becomes marginal 
in the limit of m — > 00 as y = ^-j- so that y — > 0. In this limit the fixed points are 
arbitrarily close in coupling constant space, and perturbation theory is viable in the 
whole region from the ultra-violet to the infra-red. 

We want to calculate the difference between the ultra-violet and the infra-red 
central charge Ac = cjjy — Cm. It has been argued that the infra-red conformal field 
theory of ( |4.10|) is given by the unitary minimal model M.[m — 1), as Ac in the 
perturbative limit y — > is given by Ac = -^y 3 + 0(y 4 ) [^, |J, and from ( f4.8| ) we get 
that 

dm) - dm - 1) = —r— - = - % - = ^_ + o(y 4 ). (4.11) 

v ; v ; m(m 2 -l) 2(2 -y)(A-y) 16 yy ' V ' 



7 As the secondary fields will have a conformal dimension of at least h+ 1 there will be no relevant 
secondary fields. 

8 The model ( f4lo| ) is integrable and was first studied in [Q. The only other integrable pertur- 



bations of the unitary minimal models are with (^(1,2) and 4>(2,i)- These three models corresponds 
respectively to the Korteweg-de Vries, Gibbon-Samede-Kotera and Kupersmidth equations [ jiT| . 

9 This can be seen by writing the Zamolodchikov metric in normal coordinates around the ultra- 
violet fixed point g = 0, Gij = % + 0(g 2 ), the beta-functions become in these coordinates fi l (g) = 
-V^i + 0(g 3 ). Ug^O only for and C (im>3)j = for 3 £ (1,3) then 

fij ! (g(( 1 ' 3 )) = for j =/= (1, 3) and there is no flow transverse to the (^(1,3) direction 0. 
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A general argument for all m has been given by a thermodynamic Bethe ansatz 
method in [14]] 10 . We now describe how Ac is calculated using our approximation 
method 11 , and we will compare this with the exact result Ac exac t = c(m) — c(m — 1). 

To construct the term (T(z, z)T(0, 0)) we use the Ward identities which follow 
from euclidean invariance 

d z (T(z, z)$!(a;i) • • ■ $ n (x n )) + \d z {Q{z, z)Q x {x x ) ■ ■ ■ ® n (x n )) 

n 

= 7i^2(5(z - Xi )d Xi - d z 5(z - ■ ■•$„,(!„)), (4.12) 

i=i 

here $j are primary fields with conformal dimension h{. For the correlator we are 
interested in contact terms vanish, and we get 

d- Zl d- Z2 (T(z 1 ,z 1 )T(z 2 ,z 2 )) = ^d Zl d Z2 (e(z 1 ,z 1 )e(z 2 ,z 2 )). (4.13) 

is the infinitesimal generator for scale transformation (hence its vanishing in the 
CFT), and in a renormalisable field theory must belong to the space spanned by 
the relevant and marginal fields defining the perturbation away from criticality in 
(0) 

N 

&(x) = 2nJ2P l (9)Mx), (4-14) 

i=l 

where f3 l (g) is the beta-function given in terms of the renormalised coupling constants 
g 0. ( |4.13| ) can thus be written 

d Zl d Z2 {T{ Zl ,zi)T{z 2 rz2)) = ^P\g)P(g)d Zl d Z2 (Mz uZl )^(z 2 ,z 2 )). (4.15) 

The correlator can be calculated in perturbative conformal field theory using 

the operator product expansion in the ultra-violet conformal field theory. The bare 
correlator is in the lowest order in Aq given by f3[ 12 



{<Kx)<j>(p)) M (m) + Ao / d 2 x'{<t){x)<j){ti)<f){x'))M{ »n),conn. + 0(A 



X 



2(2-3/) 



I) 

1 + Xo ^y)A(y) 2 \ 

y J 

10 In the first fixed points JA(m), m — 3,4, ...,12 were found numerically using the exact 
renormalisation group. 

11 Another way of calculating Ac follows from the proof of the c-theorem [||, where it is shown 
that R 2 ^§j = -jR 2 {TjfTjf) leading to Cardy's sum rule @ Ac = -|/ °° d{R 2 )R 2 {T£Tp. The 
correlator is normally only known perturbatively so this gives a perturbative estimate of Ac. 

12 conn. stands for the connected correlator. 
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where <fi(x) is the bare field 0(i,3)(z) and A(y) = ^\ZlfJfr(i+ y ) = 1 + 0{y 3 ). The 
operator product expansion coefficient b(y) can be calculated from a Coulomb gas 
representation of the minimal models using the formulas in 

h f,.\2 _ 16 (i-y) 4 / r(i+y/2) \ 4 / r(i- y /4) \ 3 f r(i~ y ) \ 2 / r(i+3y/4) \ _ 16 , 



3 (l-j//2)2 (l-3y/4)2 ^r(l-y/2); V r ( 1+ ^ 4 )/ V r ( 1+ ^/ K^^V/^ J 3 

Choosing the renormalisation conditions (<fi(x, g)<f)(0, g))\\ x \ =fl -i = /i 4 , the renor- 
malised correlator and the /3-function becomes 

(#(*.rt#(o, 9 )> = f i + ^^(M» - i) + o (9 



/%) = -V9 - ^b(y)g 2 A(y) + 0(g 3 ), (4.17) 

where 0(x, is the renormalised field and g is the renormalised coupling. The 
zeros of the (3— function, the renormalisation group fixed points, are thus guv — 0, 
g*m = *A(y)b(y) and therefore g e ( nA [ y ) b(y y 0) as the theory fl£7|) lies between the 
two scaling limits. From the Callan-Symanzik equation we get the running coupling 
constant || 

g(\x\) = \ij,x\ y , (4-18) 

u i- ^y )g (|z^-i) 

interpolating between g^y for |se| — > and g| K for \x\ — > oo and satisfying g([i~ l ) = g. 
Euclidean invariance allows us to write the correlator of the energy-momentum tensor 
as (T(z, z)T(0, 0)) = in terms of the dimensionless quantity R = n 2 zz. The 
differential equation ( 4.15|) then becomes 

^F(R) = ^R^mm). (4.19) 

A solution to this equation is given by 13 

_ _ ^2^2 / _ rR _ f R rR' _ \ 

F(R) = -^L. I # 2 (00) - 4i? y d# (00) + 6 J dR' J dR" (00) J + a x + « 2j R, 

where ai, «2 £ ^ The differential equation (|4.19 ) is a boundary value problem as 
F[R) is known in the scaling limits 

*M£EI:L (4 - 20) 



In the limit where i? — > the correlator (00) scales as in Ai(m) i.e. v - ., - . 
hence # 2 (00) ^ R y ^ 0, R J R dR'((j)(f)) ~ ^ -»• 0, and finally J^dR' f fl ' c£R"< 
13 These equations directly generalises to the case with more couplings. 
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R y — > for R — > 0. This sets the boundary value «i = cuv- As -F(-R) attains a 
finite value for R — > oo then all linear terms in F(R) must cancel and we set «2 = 0. 
Integrating ( |4.17| ) and inserting the boundary conditions gives 

Fi? = C[/y + - J ^ r ~ - + 27cA(y)b(y)g K - p ^ 

5 ,(3y-4 (3y-6 XN 



3(§v-i) 



Determining the function F(s) as F(s) = lij=i 9 =0( s )) the theory is fixed at the 
point of renormalisation = fi, and all the scale dependence is moved into the 
running coupling constant. The running coupling becomes with |x| 



9( S ) ~ 1 _ *A{y)b(y)g ^i _ ^ ' ( 4 ' 22 ) 

so that scale-transformations move around in the coupling constant space. The 1 
loop renormalisation group improved approximation to F(s) then becomes 

Fi{s) = c uv + —g (s^ 



2 V y- 1 

(/ , lMM /(2-ff)(3-j) (3i/-4)(3y-6) xx M ' 2 ' :>>) 

Here k = 4 as Fi(s) = F\{s 1 ^). To obtain the approximation for cir we then have 
to calculate the contour integral 

1 f e p/s 

h{p) = 7— / da—F^s). (4.24) 



2™ J Co 8 

In the limit m —>■ oo the ultra-violet and infra-red fixed points are perturbatively close 
in coupling constant space as noted above, hence F 1 (s) correctly describes F(s) in 
this limit and we should take p — > oo in ( |3.12| ) thus eliminating the contribution from 
the cut. The approximation then becomes lim p ^ 00 /i(p) = lim^oo Fi(s) = c* IR using 
( 3.15 ), and the approximation in this limit thus equals the RG improved perturbative 
result which is 

a * n 2 2 fy(2-y)(3-y) 

ACpert - c uv - c IR - — 7r{9m) 



2 — ' v y-i 



^A ( y )bWm ( + & ) (4.25) 
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Figure 5: The numerical-result cuv ~ Ii(p') against logp' for m = 14. Also plotted is 
the exact and perturbative values Ac exac t and Ac per t- The dashed line is the expected 
behaviour of 1(0) — I(p'). 



This is equal to the asymptotic form of the exact value c(m) — c(m — 1) in ( |4.11 ). We 



wish to improve this result using the approximation (|3.16|) . We rewrite the running 
coupling constant 



9{S) = i- ^) g (8 i/4-i) = TTJgW' 9 = J^Wn e ( "°°' ^ ( } 

All dependence of the renormalised coupling g are now moved into the parameter p 
setting s' = s\g\ 4 and p' = p\g\ A 

1 r p p/ s i r p p'I s ' 

h(p,9) = — / ds 6 — F^s) = — / ds'^F 1 (s') = h(p') (4.27) 

J Co s zm JC' s 

and then g(s') = j^nji- This contour integral can be evaluated for example doing a 
numerical integration or a series expansion in Fi(s'), we have done both. A numerical 



integration of Q4.27Q with m — 14 using the NAG Fortran Library is shown in figure 
|5j. In figure || the dashed line indicates the expected behaviour of the exact function 
1(0) — I(p') which is taken analogous to the curves of the free theories in figure |j. 
Writing Fi(s') as a power series 



oo 



F l (s') = ^ l (g(s')) = ^h n (sr ,A (4.28 

n=0 



and inserting this into ( 4.27 ) and using the integral representation 



1 If e 1/s ' 

l c , o ds> — W, (4-29) 



T(l + z) 2m 
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we then obtain 



An/4 



n=0 



r(l + n/4) 



(4.30) 



which is recognised as the Borel transform of order k = 4. This expression can be 
computed numerically (using e.g. Maple) by truncating to a finite n, and the minimal 
value can be found. 



m 


Ac 


Ac * 


Ac 


Ac n 
'—^^approxz 


1 1 
li 


n nnono 
u.uuyuy 


n 
U 


U.UUD4Z 


U.uuy iv 


1 9 

Iz 


n nnfioo 
u.uuoyy 


U.UUloU 


U.UUOoo 


n nn 791 
U.UU (Al 


lo 


u.uuo4y 


U.UU/45 


U.UU4o i 


U.UUOOO 


1 /I 
14 


n nn/i/in 

U.UU44U 


u.uuzoo 


U.UUoOo 


n nn/i/in 

U.UU44U 


10 


U.UUoO I 


U.UUzo i 


U.UUolU 


U.UUoO 1 


1 P. 
ID 


u.uuzy4 


U.UUZIO 


u.uuzoz 


u.uuzyo 


17 


0.00245 


0.00191 


0.00222 


0.00244 


18 


0.00206 


0.00169 


0.00190 


0.00205 


19 


0.00175 


0.00149 


0.00163 


0.00175 


20 


0.00150 


0.00131 


0.00140 


0.00150 


21 


0.00130 


0.00116 


0.00122 


0.00130 


22 


0.00113 


0.00103 


0.00109 


0.00113 


23 


0.000988 


0.000911 


0.000956 


0.000989 


24 


0.000870 


0.000811 


0.000846 


0.000871 


25 


0.000769 


0.000725 


0.000752 


0.000771 


26 


0.000684 


0.000650 


0.000671 


0.000686 


27 


0.000611 


0.000583 


0.000601 


0.000613 



In the table the obtained values of the numerical integration denoted Ac a p prox are 
listed together with the exact results Ac exact and the RG improved perturbative 
values Ac pert . For m < 11 the perturbative result will break down as Ac pert becomes 
negative and thereby violates unitarity. In [|I|] we used k = 2/y = (m + l)/2, but 
here we have chosen k = 4 as we then get a stricter bound on the exact function 
as explained in section 3.3. The results for k = 2/y and k = 4 are similar for 
all m calculated, except m = 11, 12 where they differ slightly. The improvement 
of the approximation (|3.16|) over the RG improved perturbative result is seen to be 
significant. In figure || Ac exact — Ac approx and Ac exact — Ac pert are plotted (scaled with 
m(m 2 — 1) so that all the points can be distinguished) against m. The horizontal 
axis is then the exact value. The figure shows that the approximation improves the 
perturbative results with more than a factor two. 

If we choose the maximum value k = 1 we get results which are almost identical 
to the exact values, we denote these by Ac appr0 x2, and plot in figure [7| Ac exact and 
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'Exact - Perturbative' 
Exact - Approx' 



10 12 14 1€ 



18 20 22 24 26 28 



Figure 6: (Ac ei . act - Ac approx )m(m 2 - 1) and (Ac exact - Ac per t)m(m 2 - 1) against 



m. 




Figure 7: Ac exact and Ac approX 2 against m. 



Ac approX 2 against m, the numbers are also listed in the table. In this case though 
we do not have a strict bound on the exact function I(p) as for the case with k = 4 
described in section 3.3. The exact function might still have a very small undershoot 
if for example the smallest non zero mass of a one particle state is small compared 
with p' m , i.e. Ac est ~ Ac exact . Figure |7| shows the very good correspondence between 
[ACexact and Ac approX 2, this is remarkable because the approximation is based on only 
a one loop calculation. 

5. Critical exponents for ip A in 3 dimensions 

The other example we will consider is </> 4 theory in 3 dimensions with 0(N) symme- 



22 



try 14 . We will study the Ising (N = 1) case in detail here, but the method applies 
for a general N. Quantum field theories have generic infra-red divergences in di- 
mensions lower than 4, for if 4 theory this can be seen by studying the 1PI 2 point 
function fl?||. A way of regulating these infra-red divergences is to either do an e 



expansion in e = 4 — D, or work with a massive theory. We will here do the lat- 



ter, and work in D = 3 following ]Tj|. We introduce the renormalised field, mass 
and coupling according to the conventions in |19| so the renormalised fields ip R are 
given as: cp = (Zi) l l 2 Lp R , Lp 2 = Z 2 ((f 2 ) R and m,g denotes the renormalised mass and 
coupling. The infra-red divergences are now removed and only show up in the bare 
correlators as non analytic dependence in the bare coupling |18|]. We will calculate 



two of the infra-red critical exponents v and rj, all other exponents follow from scaling 
relations 15 , e.g. Fisher's scaling relation 7 = z/(2 — rf). 
The free 2 point correlator, or propagator, is given as 



G (p,m) = / d D x <**(<p(x)<p(0)) = ~^—- 2 . (5.1) 
J p + m z 

The mass is related to the critical temperature 9 as 9 ~ m 2 , so that the free propa- 
gator can be written as 

G (p,9) = 9^h(l) (5.2) 

where h is regular at the origin and 7 = 1 and u — |. This scaling behaviour 
generalises to the interacting theory with critical exponents v and 7, hence from the 
scaling relations all critical exponents can be obtained from the propagator. We will 
now define the function F(s) above for respectively v and rj, and show that they are 
analytic in some sector with the correct scaling limits. 

Using that the massless theory for small momentum (infra-red region) is equiva- 



lent to the large momentum region at g = g c it follows from the Callan-Symanzik 
equation for the renormalised propagator that in the limit of small momentum it 
scales as G(p) ~ ^2=^7, hence defining 

F V (P) = P^I^/G R ( P ) + 2, (5.3) 



14 This model is one of the most studied models in critical phenomena and it is important because 
it shares its infra-red fixed point with a number of physical models, such as: polymers (N = 0), the 
Ising model (N — 1), super-fluid Bose-liquid (N — 2) and the Heisenbcrg ferromagnet (TV = 3). 

15 The calculation of the critical exponents in (p 4 theory are among the most precise calculations 
in quantum field theory, and they have been calculated using both the e expansion, exact renormal- 
isation group arguments, perturbative quantum field theory high temperature expansions, strong 
coupling expansions and monte Carlo simulations. These exponents describe the scaling behaviour 
of the theory in the scaling region and are universal in the sense that they are only determined by 
the infra-red fixed point and are therefore identical for all quantum field theories flowing to this 
point under the renormalisation group, irrespective of the underlying microscopic dynamics. 
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would satisfy the infra-red behaviour \im p ^ F v (p) = rj. To get the ultra-violet 
behaviour we can look at the spectral representation 

P M = _ 2 ^^iMpL +2 . )fc[f . M , (5.4) 

or simply use that G§(p) — > p 2+ m 2 i n the ultra-violet (the spectral density will only 
contribute with 5(p 2 — m 2 )), showing that lim^oo F^s) = 0, as it should be using 
the scaling relation u(2 — if) = 7. From this spectral decomposition it also directly 
follows that F v (s) is analytic for s G S(n — e') for e'<l 16 , and that 

^M-Inf 01 ^' (5-5) 
[ tor s — >• 00. 

Dimensional analysis gives us that 

dG§(p,m) dG§(p,m) ,^ n , / 9 \ . 

hence it follows that 

F v (p) = -m dG ^ m) /G*(p, m) = Z^m^-, (5.6) 

also showing that the ultra-violet value is zero as the gaussian theory is ultra-violet 
finite. We also want to compute u, it again follows from the Callan-Symanzik equa- 
tion for the n point Green function with s insertions of (p 2 that the scaling behaviour 



for low momenta at the critical point is [19 



G« s (\p,\q,m)~\°-~2^+^-i. (5.7) 

Hence in the limit p — > 



so we define 



p '— /G?i(p,q,m) = -2 + 7i--, (5. 

op V 



9Gf l (p,-p,q = 0,m) j/ iR 



FAv) = -F v (p)+A + p ^ ^ >—L/G« x fa-p,q = 0,m). (5.9) 



Using the same trick as before we see that F u (p) can be rewritten as 
(GylP. -P, 5 = 0, m)/G$(p, m)) dm 



16 Actually it will be analytic in S(tt) but to get the Borel transform F v (s) needs to be bounded 
at the origin so that F v (s) have to be bounded for all closed subsets of S which would not be the 
case with S(tt). 
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which shows that linip^oo F u (p) = 0, again because Z 2 is a constant in this limit. 
The equation ( |5.9|) shows that F u (s) is analytic in S(n — e') for all e' G (0,7r). 
This follows from writing = Gi,\Z^ Z^ an d using that G2,i(p, —p, q = 0) = 
— g^ I G2(p) and then rewriting in terms of renormalised quantities and doing the 

spectral representation. F u (s) then satisfies 

(5,d 

[0 for s — > oo. 

The functions F r) (g) and F u (g) have been calculated up to an amazing 7 loops in 
the coupling constant g. We will follow the convention of |2l| where g = -^-g and 



[3(g) = j^f3(g) giving ±1 as the first two coefficients in the (3 function; we will use 



j3,g to denote (3,g below to simplify the notation. The results are that [21 



F v (g) = 0.0109739369/ + 0.0009142223/ + 0.0017962229c/ 4 

- 0.0006536980/ + 0.0013878101/ - 0.001697694/, 
p v (g) =\g- 0.0631001372/ + 0.0452244754/ - 0.0377233459/ 



(5.12) 



3; 

+ 0.0437466494/ - 0.0589756313/ + 0.09155179/, 
(3(g) = - g + /- 0.4224965707/ + 0.3510695978/ 

- 0.3765268283/ + 0.49554751/ - 0.74968893/, 
where the beta-function has been calculated up to 6-loops ||21|| . The perturbative 



series for </ theory was in |22| shown to be Borel summable in the coupling of order 
k — 1, so we set F(s) = F(s 1 ). In [23[] it was shown using the idea from |?4j that 



the coefficients of the perturbative Green functions f(g) = Y2 n fn9 n for large n are 
given as 

f n = c(-b) n T(l + b t + n)(l + 0(l/n)) (5.13) 

where b,bi,c were calculated in p3| , p5| . b = 0.14777422 and it follows that the con- 
vergence radius of the Borel transform of f(g) is given by b s = 1/b with a singularity 
at —1/b. When b\ > we want to perform a Borel-Leroy transform of order k, which 
replaces a formal series / = J2 n f n z n with ^2 n (f n /F(l + bi + n/k)z n . Instead of ( |2.9j ) 
we then define 

r(i + M 1 fds^ /6 



this function will then have the limits lim p _ >00 /(p) = F IR and lim p ^ ^(p) = F uv as 
before. 

We introduce the scale parameter s via the exact running coupling g(s) (s —>■ 
in the ultra-violet) then we know that F ri (g(s)) and F l ,(g(s)) are asymptotic series 
in g(s). We now approximate the exact running coupling g(s) with the solution to 

dg n (s) ._ . , . 

s ds = -Pn(gn) (5.15) 
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with the boundary condition g' n (0) = 1, where (3 n {g) is the perturbative (3 function 
to the nth order. F v (g n (s)), F u (g n (s)) are then asymptotic series in s, which we will 
write as Fj(s) = Yli(Fj)i sl f° r 3 = V: v - The Borel-Leroy transform given by ( |5.14| ) 
becomes 



n+l 



As we are only working with a truncated series we will set k — 1 (only infinitesimally 
different from k = 1 + 8, 5< 1), note we are summing up to n + 1 as the nth order 
(3 function has n+l terms, and Fj are known to the (n + l)th order. 

To find the critical exponents is, rj we will here do an analytical continuation 
in the Borel transform (which we do by a conformal mapping) followed by a Pade 
approximation. The details of this calculation are given in the appendix. We get the 
following estimates 17 

v = 0.625 ± 0.004 (5.17) 

and 

7] = 0.0315 ± 0.0020. (5.18) 
These numbers should be compared with is = 0.6304±0.0013 [0, is = 0.6300±0.0015 



H], v = 0.6290±0.0025 (e expansion @), v = 0.6289±0.0008 (Monte Carlo) and for 



the other exponent r] = 0.0355±0.0025 MX, r) = 0.032±0.003 m,V = 0.0360±0.0050 



(e expansion ||), r] = 0.0374 ± 0.0014 (Monte Carlo), rj = 0.0347 ± 0.001 (strong 
coupling P7|), most of these numbers are taken from [p0[] . 

In the usual evaluation one determines the value of the critical coupling g c and 
then evaluates the re-summed series at this point. The values then becomes very 
dependent on the estimate of g c . One advantage of our method is that we do not have 
to estimate g c , likewise we do not have to perform the Laplace integral, but instead 
the critical point is reached taking the limit in the Borel transform. The disadvantage 
of this method is that we do not have specific knowledge of the quantities r c , h\ 
governing the asymptotic behaviour for the transformed series in the scale parameter. 
In the usual case of the perturbative expansion in the coupling these parameters are 
obtained from estimates of the higher order behaviour of perturbation theory. 

17 The error in v is seen to be larger than the one in 77, in contrast to the errors form other methods. 
The reason for this is that we have chosen a conservative estimate where we have averaged over the 
results with to different /3-functions (using its alternating behaviour) as explained in the appendix. 
The two /^-functions differ only in the 7th order term, and here the term is about a factor 100 
smaller than the F v term resulting in less sensitivity to this averaging. The errors obtained could 
be lowered e.g. using one of the techniques from ]2q| or as in [^9| using the information about the 
pole at r c to obtain an extra order in the Pade table, but the aim here has not been to get a very 
low error, but to apply the method described above to a well known example. 
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6. Conclusion 

We have described a general method of obtaining the infra-red limit of a physical 
quantity as an integral in the ultra-violet region. This was done using the analytic 
structure of Green functions in a complex scale parameter and by moving all scale 
dependence into the running coupling. The infra-red limit is then given as the limiting 
value of the Borel transform in the scale of the physical quantity. Changing the way 
in which the scale is introduced amounts to changing the order in the Borel transform. 

We have tested this on two examples where the infra-red limit is well known, 
namely for the central charge of the perturbed unitary minimal models and the crit- 
ical exponents of ip A theory in three dimensions. For the perturbed minimal models 
we showed how an approximation method can be obtained for calculating the central 
charge, and we get approximations close to the exact values already at one loop by 
using the largest domain of analyticity implied by spectral decomposition. For the 
(p 4 theory our estimates of the critical exponents are within the errors of other, more 
elaborate approximations. 
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A. Pade approximation 



The Borel transform is analytic in a sector S(e") where e" <C 1 jl6fl . The series given 
by ( [5.16| ) will have a pole at the negative real axis for some value p = —r c (r c > 0) 
determining its convergence radius. We will now analytically continue I(p) given by 
( |5.16| ) doing the conformal transformation (as in |T 



P 



r c + p 



P 



-oo, — r r 



P 



r c t 

1-t' 



(A.l) 



so that I(t, bi) = bi). If we write ( |5.16|) as I(p) = J2i hp 1 then I(t) is given by 

the series 



n+l i / / _ 1 \ 

1=1 m=l ^ ' 



(A.2) 



If all poles of I(p) lies in (— oo, — r c ), as argued in |26|, [HJ, then I(t) given by this series 
is convergent for t G [0, 1). We want to evaluate this expression at the convergence 
radius t = 1 (where I(t) is regular), to do this we will use Pade approximants 
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[2~8| , p9| . In Pade approximation a function is approximated by a rational function, 
called [L/M], with polynomials of degree L and M in the numerator and denominator 
(and the constant in the denominator is 1). If the n first terms of the function is 
known then these coefficients are matched with the coefficients of the polynomials 
where < L + N < n. One then forms the Pade table with entries [L/M]. The Pade 
approximation is based on the conjecture that there is a subsequence of diagonal Pade 
approximants [L/L] which converge uniformly to the function, and this conjecture 
has shown to hold in practice. The diagonal Pade approximants are conformally 
invariant and are therefore independent of r c above. Generally we should use Pade 
approximants close to the diagonal. Note that according to the Pade conjecture 
one should still expect convergence of the Pade approximants even if there are a 
finite number of poles within \t\ < 1, which is the case if not all poles of I(p) are in 
(-oo, -r c ]. 

We also have to determine the values of r c and bi. Let us first note that from 
the boundary condition g' n (0) = 1 we have that g n (s) = s + 0(s 2 ) and in Fj(s) = 
^2i(Fj)is l we thus have that (Fj)i = (Fj)i + • • • where (Fj)i is the Zth coefficient of 
g in Fj(g), this means that the convergence radius of the transformed series cannot 
be larger than the convergence radius of the series in the coupling g, i.e. < r c < 
b s = 1/6. Also we would suspect that bi ~ b' h where b\ is the Leroy parameter of the 
series in g. If bi is chosen too large we will divide by more than the actual asymptotic 



increase in ( A..2;) and the approximation value will be to small, if bi is chosen too 



small we expect to get a poor convergence in the Pade table. In the same way we 
expect the approximate value to be too small if r c is chosen smaller than the actual 
convergence radius because the values we are calculating are increasing from zero. 
This is the behaviour we see in the tables and we estimate r c and bi to be respectively 
the largest and smallest value so that there is convergence in the Pade table. 

From the asymptotic form of the (3 function ( 5.13|) it follows that it is alternating, 



and an approximation to a truncated alternating asymptotic series Ym=i fi * s to use 
the series Y^i=i fi where /; = fi for I < n and f n = f n /2 |30| . We have also obtained 



the approximation where the perturbative /3-function to the nth order (for n = 5 and 
n = 6) is approximated in this way. 

Using the criteria mentioned above, for the choice of r c , b\ we get the following 
Pade table 18 



/ ■ • 0.407663 0.405138 0.386709 0.396082 \ 

0.362383* 0.405248 0.408042* 0.392824 

0.414164 0.395111 0.392686 0.395422 

0.377688 0.392075 0.394118 

0.4843* 0.403122 
\0.406194 / 



(A.3) 



18 The numbers marked with a * has a pole close to or in the interval (0, 1). 
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for the function F u (s) leading to F u jr = 0.402 ± 0.007, where the error is the inter- 
tabular error for the points chosen (which is here larger than the Baker-Hunter error 
[j32]j ), and r c = b s , bi = 2.4. Using the ft function with 1/2 times the last coefficient 
we get the Pade table (we get the best convergence for the 5. order expression) 



/ 



V 



0.407663 0.405138 0.404239 \ 
0.362383* 0.405248 0.403736 
0.395111 0.404338 
0.400335 

J 



(A.4) 



giving F v jr = 0.404 ± 0.004, again with r c = b s and b[ = 2.4. Averaging over these 
two we get for the critical exponent v = 0.626 ± 0.003. At this value of b\ and r c we 
have the best convergence in the Pade table, a more conservative estimate of F v jr 
is obtained by varying r c and b s in a region around these values and then sample the 
highest and lowest value for which there is some convergence in the Pade table 19 . 
This gives F vIR = 0.40 ± 0.01 leading to 



v 



0.625 ±0.004. 



(A.5) 



For the other exponent rj we get approximately the same Pade tables using either 
the /3-function as given in ( 5.12Q or with a half times the last term, the convergent 
Pade table becomes 



/ • • • • • \ 

0.0291884 0.0321865 0.0326601 
0.0271003 0.0489937* 0.0327393 
0.0323937 0.0310551 0.031359 
0.0321239 0.0312882 
\0.0312289 ) 

and we obtain the critical exponent rj = 0.0319 ± 0.0010, where r c 
Again varying around these values gives the estimate 



(A.6) 



V 



0.0315 ±0.0020. 



b s and b\ = 1.8. 



(A.7) 
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